function P9 = Pding2Pzong(p4,JMV)
for i = 1:length(JMV(:,1))
    %% 构建单元压力
    pp = [p4(JMV(i,1)) p4(JMV(i,3)) p4(JMV(i,9)) p4(JMV(i,7))];
    %% 节点 1 压力
    P9(JMV(i,1)) = p4(JMV(i,1));
    %% 节点 2 压力
    kesi = 0;
    ita = -1;
    Fyp = [1/4*(1-kesi)*(1-ita);
        1/4*(1+kesi)*(1-ita);
        1/4*(1+kesi)*(1+ita);
        1/4*(1-kesi)*(1+ita)];
    P9(JMV(i,2)) = pp*Fyp;
    %% 节点 3 压力
    P9(JMV(i,3)) = p4(JMV(i,3));
    %% 节点 4 压力
    kesi = -1;
    ita = 0;
    Fyp = [1/4*(1-kesi)*(1-ita);
        1/4*(1+kesi)*(1-ita);
        1/4*(1+kesi)*(1+ita);
        1/4*(1-kesi)*(1+ita)];
    P9(JMV(i,4)) = pp*Fyp;
    %% 节点 5 压力
    kesi = 0;
    ita = 0;
    Fyp = [1/4*(1-kesi)*(1-ita);
        1/4*(1+kesi)*(1-ita);
        1/4*(1+kesi)*(1+ita);
        1/4*(1-kesi)*(1+ita)];
    P9(JMV(i,5)) = pp*Fyp;
    %% 节点 6 压力
    kesi = 1;
    ita = 0;
    Fyp = [1/4*(1-kesi)*(1-ita);
        1/4*(1+kesi)*(1-ita);
        1/4*(1+kesi)*(1+ita);
        1/4*(1-kesi)*(1+ita)];
    P9(JMV(i,6)) = pp*Fyp;
    %% 节点 7 压力
    P9(JMV(i,7)) = p4(JMV(i,7));
    %% 节点 8 压力
    kesi = 0;
    ita = 1;
    Fyp = [1/4*(1-kesi)*(1-ita);
        1/4*(1+kesi)*(1-ita);
        1/4*(1+kesi)*(1+ita);
        1/4*(1-kesi)*(1+ita)];
    P9(JMV(i,8)) = pp*Fyp;
    %% 节点 9 压力
    P9(JMV(i,9)) = p4(JMV(i,9));
end

end